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1. INTRODUCTION 

Since the discovery of the intrinsic Josephson 
effect not only the typical properties of conven- 
tional junctions were demonstrated jjj, but also 
unique and surprising features like the coupling 
of Josephson oscillations to phonons || have been 
discovered. There has also been a considerable in- 
terest in the influence of nonequilibrium effects on 
the I — U-characteristic and collective modes ||- 
^| , as the quasi two-dimensional superconducting 
layers are expected to be more sensitive to exter- 
nal perturbations than bulk materials. Recent 
measurements of Shapirosteps on the resistive 
branch |f^,|ll| suggest to study the role of non- 
equilibrium effects on these, as the in depth the- 
oretical understanding of this problem is impor- 
tant for any high-precision applications of high 
temperature superconductors (HTSC) as a volt- 
age standard. 

This paper is organized as follows: The micro- 
scopic theory for the electronic transport between 
the superconducting layers is developed and the 
close analogy of the static case to the normal state 
is pointed out. Then the consequences for the I- 
U-curve, the dispersion of collective modes and 



the position of Shapiro steps are being discussed. 
Further technical details can be found in |?],|| and 
in a forthcoming publication ]l2fl. 



2. TUNNELING THEORY 

We consider a stack of N + 1 (superconduct- 
ing) layers I — 0,1,..., N forming N intrinsic 
(Josephson) junctions in the homogeneous case. 
The normal conducting electrodes attached at the 
top and bottom of the stack in a 2-point measure- 
ment are denoted as I = — 1, N + 1 (cf. Fig. [I]). 
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Figure 1. Schematic picture of the mesa of su- 
perconducting layers in a 2-point measurement. 
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As a motivation for the following discussion let 
us first recall the situation in the normal state. 
There the electrical current (density) 

$?i = ^(^x-mD (i) 

= a u+1 Eu +1 + ^{ti +1 -ti) (2) 

between the layers I and I + 1 is given by the 
difference of the electrochemical potentials [i\ c — 
[if — cAqj (Fermi energy) in neighbouring lay- 
ers [d distance of layers, e = e|). This can 
be separated in a diffusion term driven by the 
difference A^/i^ := — [i\ of the chemi- 

cal potentials [i\ and a field term o^+iE^z+i = 

0M+i(^V; — ^o,i+i)- 

In turn, the (static) charge fluctuation 

Pi,n = -2eJV M A»i (3) 

on the layer I is determined completely by the 
filling of the conduction band or the chemical po- 
tential [il respectively. Using this and the Poisson 
equation (e^: background dielectric constant) 

Pi = eo£oo (Ei,i+i - El-l,l) (4) 
\x1 can be eliminated from equ. |^: 



„• _ ANN 
J — JLl+1 



i+ i (l - a A^A Ei <l+1 



(5) 



with the discrete derivative A^ 2 ) fi 
2fi and the coupling constant 



fi+i+fi-i- 



Qo = 



2e 2 dN 2d 



< 1 



(G) 



For a fixed dc-bias current j equ. || can be used 
to determine the electric field £7,2+1 by applying 
the operator 1 + aoA^ 2 ) on equ. M. If all con- 
ductivities 07,2+1 are equal, no charge fluctuations 
accumulate and £7,2+1 = j/a, while in the case 
of only one barrier with a higher resistance (e.g. 
00,1 <C 07 I / ) the electric field is not only 
localized at the highly resistive junctions, but is 
spread to neighboring junctions 

£0,1 = (1 - 2a )j/<7 0) i , (7) 

£-1,0 = £1,2 = aoj/00,1 . (8) 

In the superconducting state both the electric 
field £7,2+1 and the difference (/ = 0, . . . , N — 1) 

2e f l+1 

=Xl~ Xi+i 'J- J A z (z,t)dz (9) 



of the phases xi °f the superconducting order pa- 
rameters A; = |Ai|e lXl in neighboring layers are 
gauge invariant quantities and have an indepen- 
dent physical meaning. Its time derivative di- 
rectly leads to the general Josephson relation 



— 72,2+1 (i) = Vi,i+i- 



■iH+i-m z = o,.,w-i(io) 



between the phase 72,2+1 and the voltage V;,2+i = 
di t i+\Ei^ + \ across the junction, where the scalar 
potential Hi is given by: 



1 = 0,1, 



.,N 



(11) 



As can be seen from the Gorkov equation for the 
single particle Greens function, the static part 
pf c /e plays the role of the chemical potential (if 
in the superconducting state. 

In addition to this, a tunneling Hamiltonian 



H T 



^ tr c l+lra C lra e 



h Ji y ' ' +C.C.(12) 



between the (BCS-like) superconducting layers 
has to be specified explicitly in order to calculate 
jz,2+i. In the following we will assume a hopping 
matrix element f^, = t 2 s + t^g^g-^i, which has 
both a constant s and d-wave (i.e. <?k) part in or- 
der to get both a Josephson current and a finite 
quasiparticle conductivity at the same time. The 
general structure of the theory will be indepen- 
dent of this choice. 

A time dependent perturbation theory up to 
the second order in the hopping matrixelemcnt 
tkk' is performed, but no assumption for the 
external electromagnetic potentials Ao.i(t) and 
A z (t) is made. The resulting expressions for the 
tunnel current 



31,1+1 = s 



J ^ , P±(<,i')e^ (7M+l(t)±7M+l(t ' )) (13) 



between the layers / and 
fluctuation 

pi = -x P pPi + Qt 2 (t,ii,i+i) 



1 and the charge 



(14) 



on layer I are nonlinear in 72,2+1- Thereby Xpp 1S 
the charge suszeptibility of the superconducting 
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layer, which is for lu = connected with the two 
dimensional density of states N2d and indepen- 
dent of temperature T: x PP (y> = 0, T) = 2e 2 N 2 d- 
The correlation functions P± and the charge con- 
tribution Q t 2 in O(fkk') & l so depend on the scalar 
potential pi and have been expressed by Feynman 
diagrams in 

Using the general ansatz for the phase 



7m+i = 7o; + fyf + 5j lt i + i(t) 



(15) 



in resistive (f2/ ^ 0) and superconducting (f2; ^ 
0) junctions, and pi can be linearized in the 

(small) oscillations S"fi,i+i and the scalar poten- 
tial pp. 

■ ■ , ho- . 1 A(l) 

Ji,i+i = JcSin^i i+i + — — 7j j +1 - (T\-AS >fii 

2ed a 



2ed 



= -X P pHi + ~(jf!i + i-3?-i, l ) , (17) 

jg +1 = aiEi,i +1 + {ai ~ g2) ( W+ i - Mi) -(18) 

Here we restricted the discussion to a linear quasi- 
particle characteristic for simplicity. The conduc- 
tivities cTi(Cli,u) will in general depend on both 
Tifli/2e — Vdc + M/+i ~Mf c an d the oscillation fre- 
quency lu of Sjij+x, which coincides with fl e.g. 
on the autonomous first resistive branch. Prod- 
ucts like (<r E l>l+1 )(jt) = / a (t - t')E u+1 {t')dt' 
are to be interpreted as a folding in time space. 

All conductivities <Ji are of the same order of 
magnitude as the experimental quasiparticle con- 
ductivity er cxp w 2(fc£!cm)~ 1 (Bi-2212) and con- 
sequently corrections <Ti/e Lu <C 10 -2 are small 
for frequencies to > to p \ on the resistive branch. 
The conductivity a\ is presented for different val- 
ues of fl in fig. g, There the transition to the 
conductivity <jn — 4e 2 7rA f 2 dis/^ 3 f° r w > 2eA 
(Ao: maximal d-wave gap) and the slight differ- 
ence between the value of o\ in a superconducting 
and a resistive junction can be seen. The nega- 
tive part for small frequencies is not a violation of 
causality, but arises from the fact that a\ is not 
a conventional transport coefficient. 



The static limit of <j\ is subtle and deserves 
special attention. Formally the result for lu = 
cannot be obtained in our formalism by expand- 
ing pi or in the static part pf° of the scalar 
potential directly. Instead of this the static in- 
homogenity has to be included in the equilibrium 
Greens function in zeroth order in tkk' ■ For lu — > 
and T < eV <C A one gets the finite value 



<7l(f2, LU — > 0) ~ (TJV 



T eV 
A^A^ 



(19) 



while the strictly static function af c := ct\{lu = 
0,T) = vanishes exactly for all temperatures. 
This discontinuity is due to the fact that our 
only relaxation mechanism is the incoherence of 
the hopping £kk' between layers. Therefore in 
a tunneling process (which is already of order 
^kk') 110 additional relaxation within the layers 
is possible, if we work strictly in 0(t^ k ,). In- 
cluding strong impurity scattering in the layers, 
which is the main mechanism for charge imbal- 
ance relaxation in <i-wave superconductors, leads 
to <7\{lu — * 0) = <J\(u> = 0) = like in ||, which 
will be considered mainly in the following. It is 
pointed out that this choice is strictly speaking 
model dependent and a finite value af c similar 
to our result for o\ {to — > 0) in cqu. [l^ could be 
obtained, if additional scattering mechanisms are 
considered. 

The elimination of V/^+i, pi and Spi from equa- 
tions |, |TJ, [TJ and U leads to a set of coupled 
equations of motion for the phases 7;,z+i (lu^ := 

2edj c /he , lu c = u^eq/ctq, A^ 2) := 1 - A^a): 

i = A™ sin 7M+1 +— A( 2 ) 7u+1 + ^Af } 7 U+1 (20) 

The coupling constants are given as: 



a(u>) 



dXpp 



1 



CM := - 



e €LU 

i(J2 2<7i 

e e «cr 2 



(21) 
(22) 
(23) 



These are in the static case completely deter- 
mined by the normal state value in equ. ||: 
a(tu = 0) = rj(u = 0) = ao and £(u> = 0) = 0. 
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Figure 2. Real and imaginary part of o\ in 
superconducting (f2 = 0, top figure) or resistive 
(ft ^ 0, bottom figure) junctions at T = 0. 



The rough theoretical estimate for an ~ 0.28 
(for e = 10, d = 1.5nm) based on a twodimen- 
sional electron gas with density of states A^d = 
m e i/2-!T 2 h 2 (not including spin) is to be considered 
as an upper bound, as the density of states at the 
Fermi surface could be enhanced by a factor 2-3 
due to strong electronic correlations J]l| . A more 
reliable estimate a « 0.2 is possible from opti- 
cal experiments and will be elaborated in detail 
elsewhere 

Linearizing the set of equations ^0] one ob- 
tains the dispersion of collective modes (in order 
0(*kk')> := CT /e e = w^/u c , uJ 2 pl = w^/ef ): 



2 H5 1 -2 



1 

pi r 



1 + i]k 2 z 



<*? -»ttM= (24) 

57^ (l + a/c^) - ioJojQ (l + ?jfe^) (.25) 

For frequencies cu p \ <C lu <C Ao on the resistive 
branch these represent weakly damped plasma 
modes of the superconducting condensate. For 
the special case e — > the dispersion u>(k z ) —> 
for k z — > of the Goldstone mode associated with 
the spontaneous breaking of the gauge symmetry 
in the superconducting state can be obtained ex- 
plicitly. This reproduces in the Kuboformalism 
the results in which were obtained by solv- 
ing kinetic equations. 

Taking into account only the leading order 
0(i£. k , ) in the hopping matrix element, a(w,T) 
has been calculated in 0] and is presented for 
T = in fig. ||. This result is obtained for a 
constant background dielectric constant e, while 
the frequency dependence of 600(0;) near phonon 
resonances can change the behaviour of a(uj) con- 
siderably (^] . For the temperature dependence see 
§• 

We would like to point out the remarkable sim- 
ilarity of the static quasiparticle current in the 
superconducting and in the normal state. Due to 
u\ (id = 0) = cqu. |l6| reduces to: 



Jdc = 3c{sin^i +1 )+—{^i M 



,0-1 



}+(-±( w+1 _ w ))(26) 



In the last term of equ. 
does not contribute due 



26 the static part of \x 
ax(Q,u = 0) 



tO (Ti(ll,UJ = 0) = 0, 

but a finite dc-contribution might arise from the 
combination of /i(^) and Josephson-likc terms in 
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Figure 3. Real and imaginary part of a(u>) for 
d-w&vc pairing for temperature T = in O(^k') 
(lu in units Aq(T), — const.). 



o\. This term is in intrinsic junctions negligible 
of order at^, -C 1 and vanishes exactly for SN- 
junctions: 



Jdc 



j c (sin7 M+ i) + ^(<y M+x ) = 
j c (sm lu+1 )+a a A^E^ +1 



(27) 

(28) 
(29) 



where the quasiparticle current coincides with 
equ. 0, if we express "fi,i+i by £^+1 and /x; via 
equ. M The dc-equations 



J 



3 



rr SN rr 



SN 



SN\ 



a\- ) A (1 Vjv = 
- a' E N - hN ) , (30) 
(31) 



r SN 



SN 



((1 + oi' )En,n+i 
{(l + a' )E_ lfi -a' E ,i) 



total voltage V dc = J2i=-i Vu+i ^ s exactly given 
like in a stack of independent junctions, i.e. the 
coupling constant ao does not enter in the I- 
V-curve. The dc-component of the supercurrent 
might add a correction to this due to the interac- 
tion a at the frequency f2, but this contribution is 
suppressed by a factor t^ w ,ao <C 1 and therefore 
no significant effect on dc-properties are expected. 

It is pointed out that this general result cannot 
be obtained by assuming an ohmic quasiparticle 
current j ; q f +1 = aEij+i as in An additional 
indication that the correct choice of the quasi- 
particle current is essential will come from the 
discussion of Shapiro steps in the next section. 

This might also be the reason that up to now 
no clear experimental evidence for the coupling a 
in the T-V^-characteristic, like small deviations in 
the additivity of resistive branches ||, could be 
found. 

Nevertheless, the distribution of the electric 
fields Efi +1 within the stack is affected by the 
coupling ao like in equ. and || in the normal 
state. For a single resistive junction the elec- 
tric field is not contained in this junction, but 
leaks into neighboring junctions in the same way 
as around a junction with much lower resistance 
in a stack of normal junctions. 

It is stressed that the above discussion makes 
use of the fact that af c — in the presence of 
strong impurity scattering in the layers. In a 
more general model this might not be the case 
and coupling terms ~ aoaf c /o~Q C appear. 

The coupling via charge fluctuations might also 
have a considerable effect at frequencies larger 
than the charge imbalance relaxation rates (e.g. 
for phaselocking at microwave frequencies), where 
the cri/co in equ. ^ can be of order O(l) in any 
model. 



for the SN-junctions are slightly different, as 
the change of the chemical potential /j, c l 1 = 
l-i c N+ i = in bulk materials is negligible due 
to the large density of states. Thereby a' = 



a (l 



'l,SN 



/°0 sjv)j wri ich is in the model with 



strong impurity scattering in the layers given as: 



a . 



If we neglect the supercurrent for a moment, 

(2) 

the inversion of the operator Aq q shows that the 



3. SHAPIRO STEPS 

A Shapiro step is generated by microwave ir- 
radiation of frequency f2, which creates an os- 
cillating component of the c-axis bias current 
-fac = Iii sin(fit) . Thereby the dynamics of the 
phase in a resistive junction is given as 



Jn,n + l(t) = 70n + &t + 7i„ SUl(£l£) 



(32) 



G 



Therefore the supercurrent (J^: Bessel functions 
of first kind) 



J c sin(7 n . n+ i(t)) = 

Jfc(7ni)sin(7„ - 



E 

k= — oo 



(i 



has a dc component 



J? 



I c J-l(7l„)sin(7 0n ) 



(33) 

k)nt) 

(34) 



which creates the (first) Shapiro step in the I- V- 
curve. 

The nontrivial question in the following will be, 
what dc-voltage V^tep is connected with a given 
f2, if the generalized Josephson relation equ. 10 
is taken into account. Using equation |] and 17 
this can be presented as: 

= (l+2a)E u+1 ~a(E^ u + E l+1J+2 )(35) 

Zed 

for I — 0, . . . , N — 1. This set of N equations 
is not sufficient to determine all N + 2 voltages 
VJ^+i and has to be complemented by the current 
equations ^ and [3l] in the SN-junctions. 

Note that for the total voltage V s t cp across the 
stack the transport coefficients in the intrinsic 
Josephson junctions are irrelevant. 

In the "superconducting" state, which is de- 



fined by (7,2+1} = for all I 
easily obtains (in O(a' )) 



0,. 



, N — 1 , one 



1.0 



e: 



dc 

1,1 + 1 

.dc 
0.1 



E\ 



= E 
= 
= E 



N,N+1 



(1 



a 



1 = 1, 



dc 

N-l.N 



0/j/°0,SN 

N -2 , 

dc 



(36) 
(37) 
(38) 



and the total dc-voltagc 

V dC = E <^ = ^ 09) 

l = -l a 0,SN 

This linear branch reflects the contact resistances 
at both SN-junctions and is not affected by the 
coupling a' a , although the voltages V^+i in the 
intrinsic junctions in general do not vanish. 

For a single resistive junction in the top layer 
I = N - 1 ((tv-i,jv) = ^, (7n,n+i) = 0, n else) 
one obtains analogously: 



E 



N,N+1 



(l-<)j/<SN + <^ d + 0(a'i), 



stop 




T dc 
T 1,SN 



i-^o-i— = (40) 



_dc 



(41) 



The last result is correct in all orders of a' . Here 
the dc-conductivity &f c SN = has been kept ex- 
plicitly in equ. ^ in order to demonstrate that 
the position of the step in general depends on 
the value of crf c SN . In the model considered here 
crf c SN = and the Shapiro step is where ex- 
pected, which would not be obtained by using 
the ohmic quasiparticle current alone (cq c sw = 
<T iSAf) like in H. Nevertheless this result opens 
up the principal possibility that there is a shift in 
the position of the Shapirostep in a more general 
theory (e.g including electron-phonon scattering) 
due to nonequilibrium effects, where af c might 
not vanish. Taking out result equ. [19 seriously 

fl-5 Km' 



the relative shift would be small 
tectable. 
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Figure 4. Schematic 4-point geometry by |lf 



This result is also valid in a 4-point measure- 
ment geometry, as sketched in Fig. |] and used by 
PI . Although the total voltage 



^stcp 



N ~ 1 fcO 



1=0 



\ . 2ja a 
a o) + Tdc— 



(42) 



'0,SN 
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across the intrinsic contacts depends on ao, this 
deviation of the expected position of the Shapiro 
step cannot be detected at the contacts. The elec- 
trochemical potential is constant along the super- 
conducting layers, but the potential fJ,i(x) and the 
electric field E\\ are not and will compensate the 
additional contribution in equ. [f|]. Also in the 
non equilibrium theory of || there is no shift of 
Shapiro steps. 

Finally, note that all the above considerations 
only apply to an experimental situation where 
the electric field E dc is measured directly to de- 
termine the dc- voltage V dc . This might not be 
operational with conventional voltmeters, which 
actually detect a current through a circuit with 
high resistance, which is driven by the difference 
of the electrochemical potentials fi° c rather than 
the voltage V dc @. 

4. CONCLUSIONS 

The microscopic theory of tunneling in layered 
superconductors has been studied including the 
effect of the scalar potential on the layers and an 
estimate of the coupling constant a « 0.2 has 
been given from optical experiments. The static 
quasiparticle current includes both field and dif- 
fusion terms as in the normal state, which turns 
out to be crucial for all dc properties. In the 
model with strong impurity scattering in the lay- 
ers, Shapiro steps are exactly at the expected po- 
sition T4tcp = M2/2e. In this case the effect of 
charging on the total I-V-curve is suppressed by 
a factor ao^kk' an( ^ therefore negligible, but it af- 
fects the distribution of the electric field within 
the stack. 

This results motivate the study of more general 
microscopic models including different relaxation 
mechanisms like electron-phonon scattering, as a 
finite af c would modify the position V stcp in a 
characteristic way with important consequences 
for potential applications as a voltage standard. 
Also precision experiments on Shapiro steps could 
then be used as a sharp test for microscopic the- 
ories. 
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